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Abstract 

We present an exact solution of the 0(n) model on a random lattice. The coupling 
constant space of our model is parametrized in terms of a set of moment variables and 
the same type of universality with respect to the potential as observed for the one- 
matrix model is found. In addition we find a large degree of universality with respect 
to n; namely for n g] — 2, 2[ the solution can be presented in a form which is valid not 
only for any potential, but also for any n (not necessarily rational). The cases n = ±2 
are treated separately. We give explicit expressions for the genus zero contribution to 
the one- and two-loop correlators as well as for the genus one contribution to the one- 
loop correlator and the free energy. It is shown how one can obtain from these results 
any multi-loop correlator and the free energy to any genus and the structure of the 
higher genera contributions is described. Furthermore we describe how the calculation 
of the higher genera contributions can be pursued in the scaling limit. 



1 On leave of absence from NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 



1 



1 Introduction 



The 0(n) model on a random lattice JTJ, [|] is a matrix model which regarding its 
complexity can be placed somewhere in between the one-matrix model and the two- 
matrix model. It is therefore a natural intermediate step if one wants to study the 
generalization of 1-matrix model techniques and results to the two- and eventually the 
multi-matrix case. The model is also interesting in its own right having an appealing 
geometrical interpretation and a very rich phase structure || [!§, f|, [|. In particular 
when n = 2cos(utc) with u — l/k, < I < k and l,k £ Z the model has critical 
points for which the associated scaling behaviour is that characteristic of 2D gravity 
interacting with rational conformal matter fields of the type (p, q) = (k, (2m +l)k±l) 
and with v general any central charge between c = — oo [v = 1) and c = 1 (v = 0) can 
be reached. However, the continuum theories that one obtains from the 0(n) model in 
the rational case contain only a subset of the operators of the corresponding minimal 
models. |§ §. 

In the present paper we will solve the model exactly, i.e. without any assumption of 
being close to a critical point. The genus zero contribution to the 1-loop correlator will 
be calculated solving the saddle point equation of the model, following the idea of refer- 
ences || and and the higher genera contributions by a generalization of the moment 
technique of reference [0. As usual this technique will allow us to find from the 1-loop 
correlator any multi-loop correlator as well as the free energy. The parametrization of 
the coupling constant space of the model in terms of moment variables reveals that 
the model possesses the same kind of universality with respect to the potential as the 
one-matrix model. In addition there appears a large degree of universality with respect 
to n. 

In the case of the one-matrix model the moment description facilitated the analysis 
of the double scaling limit M. For example the result that the continuum 1-matrix 
model partition function is a r-function of the kdV hierarchy || could easily be under- 
stood in this description ||10|| , the analysis relying on a representation of the r-function 
as a matrix model, namely the Kontsevich model [11], |I2"|| , and the moment description 
of this model [13|]. The r-functions of the kdV p hierarchies with p > 2 can also be 
represented as matrix models, namely as generalized Kontsevich models and recently 
the appropriate moment description of these models has been found |H , HJ. Hence it 
should be possible to determine which is the precise relation between the continuum 
partition function of the 0(n) model, for n rational, and the r-functions of the kdV p 
hierarchies by comparing the moment description of the 0(n) model with the moment 
description of the generalized Kontsevich models. This requires of course that a d.s.l. 
relevant version of the moment description is developed for the 0(n) model. A part of 
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our paper will be devoted to the development of such a description. 

We will start by, in section 2, presenting the model and the most important equa- 
tions needed for its solution. Then we will proceed with the exact solution, for 
n g] — 2,2 [ in section 3, and for n = ±2 in section 4. Section 5 is devoted to the 
study of the double scaling limit and section 6 contains our conclusion and a discussion 
of possible future directions of investigation. 



2 The Model 

In the following we will consider the O(n) model on a random lattice, given by the 
partition function 



Z = e N2p = [ dM TT dAi exp ( -N Tr 

JNxN H[ V 



n 

2 



8=1 



(2.1) 



where M and A i} i = 1, . . . ,n are hermitian N x N matrices and 

V{M) = f] —M j . (2.2) 
j= i 3 

In the language of Feynman diagrams the model describes a gas of n different types of 
self-avoiding loops; non- interacting and living on a random surface?] 0. To begin with 
n is an integer but by analytical continuation the model can be defined also for non 
integer values of n. We will restrict ourselves to the case |n| < 2 and we will use the 
following parametrization 

n = 2cos(z/?r), < v < 1. (2.3) 

We note that for n = the model is identical to the usual hermitian 1-matrix model. 
Furthermore for n = 1 and a special cubic potential the model describes the Ising model 
on a random lattice || . We shall in particular be concerned with the calculation of the 
free energy, F, and correlators of the M-field of the following type 

W( Pl , . . . , P s) = N s - 2 ( Tr Tr . . . Tr -^-tt) (2.4) 

\ pi-M p 2 -M Ps -M/ conn 

The genus expansion of these objects reads 

F = £ N- 2 *F 91 W( Pl , . . . ,p.) = ^ N-^W'fa, ..., Ps ). (2.5) 

9=0 9=0 



2 Strictly speaking, to have this interpretation, we should include mass terms for the A-fields and 
exclude the term linear in M in our action. However, this rearrangement can be obtained by a linear 
shift of the matrix M and since we will work with a generic potential such a shift can always be 
performed in the final result. 
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and we have 

W '^-^ = l^-W^) F " S>lors>2 (2.6) 

where 

dV(p) p*+i d 9j ■ y ■ 1 

In the remaining part of this section we shall introduce the tools which will allow us 
to determine, for any potential V(M) and any n G [—2, 2], W 9 (p\, . . . ,p s ) and F g for 
(in principle) any g and any s. Eventually it will be convenient to treat separately the 
cases n g] — 2, 2[ and n = ±2 but here we shall address the aspects which are common 
to all values of n. 

2.1 The saddle point equation 

The integration over the A matrices in our partition function ( [2.1|) is gaussian and can 
directly be carried out. This leads to the appearance of a 1-matrix integral in which we 
can diagonalize the matrices and integrate out the angular degrees of freedom. By this 
procedure our partition function (up to a constant) turns into the following integral 
over the eigenvalues of the matrix M 

zoc/ 00 n^v- jv ^ y(Ai) n(^+^)~ n/2 n(^-^) 2 - ( 2 -8) 

In the limit N — > 00 the eigenvalue configuration is determined by the saddle point of 
the integral above 0. The corresponding saddle point equation reads 

As usual this discrete equation can be transformed into a continuous one by introducing 
corresponding to the matrix M an eigenvalue density p(X) = jrJ2i^(^ ~ which in 
the limit N —>■ 00 becomes a continuous function |J. When one of the eigenvalues 
approaches the origin, the integral (|2.1|) ceases to exist (cf. to equation (|2.9| )). Therefore 
we will always assume that the eigenvalues are confined to the positive real axis. More 
precisely we will consider the situation where the eigenvalue density has support only 
on one interval [a, b] on the positive real axis and is normalized to one, i.e. 

supp p(A) = [a, b], a > 0, (2.10) 

b 

p(X)d\ = 1. (2.11) 
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Of course the results obtained in this situation will allow an anlysis of the case a — > 0. 
In terms of the eigenvalue distibution the saddle point equation ( |2.9| ) reads || 

HA) = 2/V^-n/V^- (2.12) 

J a A — jl J a A + /! 

The saddle point equation can also be written in terms of the genus zero one-loop 
correlator ||. One has 

W°(p) = (2-13) 

J a p — jl 



and the conditions ( |2.10| ) and ( |2.11| ) on p(A) are equivalent to demanding that W(p) 



is analytic in the complex plane except from a cut [a, b] and that 



W(p)^-, p^oo. (2.14) 
p 



The inverse relation to fl2.13|) reads 

p(A) = — \W°{X-iO) -W°(\ + i0)\ (2.15) 

and the saddle point equation for p(A) turns into the following equation for the genus 
zero contribution the one-loop correlator 

V'(p) = W°(p + i0) + W°(p-iO) +nW°(-p), p e [a,b]. (2.16) 
2.2 The loop equations 

The loop equations of the model can be derived in various ways [§]]. Here let us use a 
formulation which exposes very clearly the analogy with the 1-matrix model case. First 
we exploit the invariance of the partition function ( |2.1| ) under the following redefinition 
of the field M 

M^M + e — - — . (2.17) 
p-M v ' 

Introducing this shift in Q2.1|) gives rise to the following equation 

i 'W^ w ^ ] + nx(p) = (w(p)) " + w>ww w m < 2 - 18 > 

where 

x(p) = ^(Tr^4 2 ) (2.19) 

and where the contour G\ encloses the cut [a, b] of W(uS) but not the point u = p. We 
will use the convention that all contours are oriented counterclockwise. Next, let us 
consider the following redefinition of the field Ai 

Ai^Ai + e — — Ai — — . (2.20) 

p-M -p-M K J 
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Inserted into ( |2.1p this shift leads to the following identity 

- x(p) - X(~P) = W(p)W(-p) + ±^—W{-p). (2.21) 



From ( |2.18p and ( J2.21 ) we can obtain a closed equation for the 1-loop correlator of the 



M-field 

= (W(p)) 2 + (W(-p)) 2 + nW(p)W(-p) 

+ {w^ WM + wh) w( - p) + n wn) w( - p) } ■ (2 ' 22) 

This equation exhibits a strong similarity with the equation for the 1-loop correlator 
of the hermitian 1-matrix model but as opposed to the latter it is non local. However, 
as we shall see in section |3J], [4.1.1| and [4.2.1| there exists an efficient way to deal with 



this non-locality. 

3 The case n e] - 2,2[ 
3.1 Reformulation of the loop equation 

With the aim of reformulating ( |2.22|) as a local equation, let us introduce a function 
W r (p) by 

= 2^(ri-»r(-p) 
4 — n z 

Furthermore, let us corresponding to a function or an operator h(p) define h±(p) by 
, . , e + ^/ 2 h(p) + e-^h(-p) 

h+ip) = — \ , h -(P) = h+(-p). 3.2 

2 sm(z/7r) 

Inversely we then have 

h(p) = -i (e iu * /2 h + (p) - e- iU7T/2 h_ (p)) . (3.3) 
Introducing the transformation ( |3.2| ) into the loop equation ( |2.22| ) one gets 

<L "' " {W r+ (u)W-(u) + W r -(u)W+(u)} 



c 2 2iri p 2 — to 2 

= w + (p)w.(p) + w 2 d^d) w - ip) (3 ' 4) 

where the contour, C^, now encircles [a, b] as well as [—b, —a] but not the point u = p 
and where d/dV+(p) is shorthand notation for (d/dV(p)) + . Introducing two linear 
operators K + and K_ by 
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and inserting the genus expansion of the correlators in Q3.4| ) we find 
{K + - Wl(p)} Wi{p) + {k- - W°(p)} wiip) 

= £ wi(p)wr 9 \ P ) + -J—wr 1 ®, g > 1. (3.6) 

The similarity with the corresponding equation of the hermitian 1-matrix model ap- 
pearing in reference [[jj is striking and we will later show how the strategy of reference 0] 
for solving the loop equation genus by genus can be generalized to the present case. 
Of course an iterative procedure for solving (|3.6|) requires the knowledge of W°(p). 
In the next section we will show how one can write down a closed expression for this 
correlator, i.e. an expression which is valid for any potential V(M) and any n e] — 2, 2[. 

3.2 W°(p) in terms of an auxiliary function G(p) 

To determine the 1-loop correlator at genus zero we follow the idea of reference || . As 
mentioned earlier we restrict ourselves to the one-cut situation. Our starting point will 
be the saddle point equation ( [2.16| ) which together with the boundary condition ( |2.14| ) 



determines uniquely W°{p). Let us split W(p) in a regular part W r (p) and a singular 
part W s {p) 

W(p) = W r (p)-W s (p). (3.7) 
From ( |2.16 ) it follows that W r (p) is given by ( |3.1| ) while W®(p) obeys the homogeneous 



saddle point equation and the boundary equation 

W°(p) ~ W r (p) - -, p^oo. (3.8) 

In the language of the rotated functions W±(p) we have the following situation 

Wl{p) = W T± (p)-W° ± {p) (3.9) 

with W r+ {p) being given by 



W r+ {p) = i ^ (3.10) 

4 — n z 



and with W®±(p) obeying the equations 



The boundary equation (|3.8|) translates into 



W° s± {p - iO) = -e ±iun W° (p + iO). (3.11) 



WSM ~ W rM - - (3.12) 

In order to obtain a closed expression for W®±(p) we introduce an auxiliary function 
G{p) with the following properties 
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1. G±(p) fulfill the equations ( J3.ll ) . 



2. G(p) is analytic in the complex plane except from a cut [a, b] and behaves as 
(p — a)~ l l 2 (p — in the vicinity of a and b. 

3. G±(p)~±j, p^oo. 

In section |3.3| we will show that these requirements are enough to fix G(p) uniquely. 
For the moment let us note that 1-3 imply that the function R(p) = G+(p)G-(p) is 
even, behaves as 1/p 2 as p — > oo and can have no singularities except for single poles 
for p = ±o, ±6, i.e. 

R(p) = G + (P)G-(P) = p l P ^l b 2y ( 3 - 13 ) 

We will choose the convention that +e is a root of G + (p) while — e is a root of G-(p). 
We will later write down an equation which determines e in terms of a and b. Now if 
we write a generic solution of ( |3.11|) , S±, as 



S±(p)=S ± (p)G ± (p) (3.14) 

we have 

S± {p - iO) = 5 T (p + zO) (3.15) 

which means that the even function S + (p) + S-(p) is a regular function while the odd 
function S + (p) — S^(p) has a square root branch cut [a, b}. Hence we have 



S + (p) = A(p 2 )+pB(p 2 )^, VP = vV - fl2 )0° 2 - &2 ) ( 3 - 16 ) 

with A(p 2 ) and B(p 2 ) regular but not necessarily entire functions. Since +e is a root 
of G + (p), A(p 2 ) and B(p 2 ) may have a pole for p = e without S + (p) becoming singular 
there provided the accompanying pole for p = — e is cancelled (cf. equations ( |3.13| ) 
and ( |3.14| )). Since we are not particularly interested in solutions which vanish for 
p = e a more convenient parametrization is 

S+(p) = A(p 2 )G + (p) +pB(p 2 )g + (p)G + (p) (3.17) 

where 

= 4^ (3-18) 
p z — e z 

and where A and B are again regular but not necessarily entire functions. We can also 
write 

S(p) = Aip 2 )Gip)+pBip 2 )Gip) (3.19) 
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where 

G(p) = - (e^ 2 g + (p)G + (p) + e -^V(p)G-(p)) (3.20) 
We draw the attention of the reader to the equation ( 3. 19]) . This equation will play a key 



role throughout the paper. It says that any solution of the saddle point equation ( 3.11|) 
can be parametrized in terms of the two functions G(p) and G(p). In particular one 
has that any such solution can be parametrized in terms of any two other independent 
solutions. A study of the analyticity properties of G(p) reveals an interesting symmetry 
of the model. Let us for a moment write the function G(p) as G y (p) where the index 
v is the parameter which enters the relation n = 2cos(z/7r). Then we have 

G"(p) = G {1 - U \p). (3.21) 



This follows from the fact that G(p) is a solution of the saddle point equation ( 3.11|) 



with v being replaced by 1 — v. Furthermore from ( |3.21|) it follows that the parameter 



e = e„ entering the relation ( [3.13|) for G u {p) is related to the corresponding parameter, 
e x - v , for G l ~ u (p) by 3 

ab 

ei-u = 3.22 

e v 

Let us now specialize to the 1-loop correlator. Since we want this function to be 
finite in the limits p — > a, b we choose in this case a slightly different but equivalent 
parametrization, namely 

W° s+ (p) = Vp(A(p 2 ) 9+ (p)G + (p) +pB(p 2 )G + (p)) . (3.23) 

Due to the assumptions concerning the analyticity properties of the 1-loop correlator 
A(p 2 ) and B(p 2 ) must here be entire functions and from the boundary condition ( |3.12|) 
it follows that they are necessarily polynomials. Using the relation (|3.13|) one easily 



concludes that A(p 2 ) and B(p 2 ) can be expressed in the following way 

A{p 2 ) = ±(G„(p)W° + (p) + G + (p)Wl(p)), (3.24) 
pB(p 2 ) = Ug-(p)G-(p)W! + (p)-g + (p)G + (pWl(p)). (3.25) 



The fact that A{p 2 ) and B{p 2 ) are polynomials and that W° + (p) ~ W r+ (p) + 0(1/ p) 
allows one to conclude 

A(p 2 ) = Even polynomial part of W r -(p)G +(p) 

^uW r ^)_ G+{u) (3 26) 



2%i p 2 — tu 2 



3 The choice of sign in this relation is a rather technical point. It relies on the expression ( |3.34| ) in 
the next section. We note, however, that all physical quantities depend only on e 2 . 
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-P f —^Fri9 + {p)G + {u) (3.27) 



pB(p ) = —Odd polynomial part of W r -(p)g + (p)G + (p) 
du W r -(uj) 
loo 2ni p 2 — oS 1 

where § ao means integration along a contour which encircles oo. In total one can write 
W®(p) as the following contour integral 

W%{p) = V^G+ip) 1 ^f^G + (u) {ug + (p) - pg + (u)} (3.28) 

where the contour C 2 encircles the cuts [a, b] and [—b, —a] of G + (u) but not the point 
to — p. The points a and b are determined by the following two equations 

/ ^( W )</ + (w)G + (w) = ) (3.29) 
Jc 2 Awi 

— W r Juj)G + (uj)uj = 7 1 ; x (3-30) 

c 2 2vri r y J +y J 2cos(z/vr/2) v ; 



which follow from the boundary condition ( |2.14p . We note that by using the analyt 



icity properties of the various functions entering the integrand ( |3.28| ) one can replace 
§ C2 tWr^uj) {. . .} by -Jggj § Cl ^V'(uj) {. . .} in the expressions (p§-(pg). It 



is a matter of taste which expression one prefers to work with. The former reflects 
more clearly the structure of the loop equation while the latter expression ressembles 
more the one of the hermitian 1-matrix model. 

3.3 Determination of the auxiliary function 
3.3.1 General case 

One can derive a differential equation for G + {p). To do so one first observes that the 
function (^y/p~G + (p)j will fulfill the equation ( |3.11|) and hence have a parametrization 
of the type (|3.17|) . When supplemented by the boundary condition for G + {p) this 
observation allows one to conclude 

^ (VpGM) = f a - ^ +pg + (pYj G + {p) (3.31) 

where a is some yet not determined constant which has the following role 

G+ ( P )= l -{l-^ + o(^y p^oc. (3.32) 

For given e and a the equation (|3.31|) determines G + (p) uniquely. It is easy to see that 
G + (p) given by the following elliptic integral is the unique solution we seek 

\Jp 2 — e 2 
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provided the following two equations hold 





i{l- V )* r dx ( e^e +a ] = ™ (3 . 34) 



These equations ensure that G + (p) has the correct asymptotic behaviour as p — > oo and 
that G + (a + iO) = -e lU7T/2 G^(a-iO) (cf. to equation (|3~TT1) ). Together they determine 
the unknowns e and a. In particular it can be shown that e must necessarily lie on the 
positive imaginary axis and behave as a? when a — > 0. One can derive another set of 
equations which determines these two quantities and which will be of importance for 
the analysis in the following sections. Using the same strategy as for the derivation 
of ( |3.31| ) one finds the following expression for the derivative of G + (p) with respect to 



a 2 



d^ G M = -rzr^ ( x aP9 + (p) + 5(1 - Pa)) g+(p) (3.35) 

where 

1 de 2 e 2 - a 2 a 2 de 2 

Xa = -2w^r> pa = (3 - 36) 

Now comparing the expressions for ^2-§^G + (p) and ^-^G + (p) that one obtains 
from ( p.31|) and ( p.35|) respectively, one finds the following relation between a, b, e 
and a 

ey^ , de 2 a 2 (b 2 - a 2 ) 

a = -^ + ^^^' (3 ' 37) 



da 



-X a - (3.38) 



da 2 

In particular these two equations allow one to write down a second order differential 
equation for e(a,b). We shall refrain from doing so since we have not been able to 
extract any further information about the model from the resulting equation. Let us 
for later convenience note that we have also the relation 

d_ 

dp 

as well as 



(Vp~9 + (p)g + (p)) = (<yg + (p)+p) g+(p) (3.39) 



() (pG+(p)) = -P^>4~, (9+(p)G+(p)) + X a g+(p)G + (p) (3.40) 



da 2 e 2 — b 2 da 2 

Here ( |3.39|) follows immediately from (|3.21|) by noting that for g + (p)G + (p) the param- 
eter a entering ( p.32|) is replaced by a, 

-a = a-^ (3.41) 

e 
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and to derive the relation ( 3.4Up one makes use of the fact that any two solutions of 



the saddle point equation can be parametrized in terms of any two other independent 
solutions. The detailed nature of the parametrization follows from the analyticity 
properties and the asymptotic behaviour of the functions involved. Needles to say that 
relations similar to ( |3.35| ) and ( |3.40| ) concerning the differentiation with respect to b 2 
follow from these by the interchangements a <-> b and that e and a depend on a and b 
in a symmetrical manner. As we will show in the next section when v is rational G+(p) 
can be further explicited. 

3.3.2 Rational case 

Let us parametrize v in the following way 

17 = J 0</<g, l,qeZ + (3.42) 
and let us following reference [|| introduce the function 

Tip) = \ {(G + (p)) q + (-l) q+l (G4p)) q } • (3-43) 

From the requirements 1-3 on G(p) it follows that T(p) is a rational function with 
poles at ±a and ±6 of order [q/2] (the integer part of q/2). Furthermore from (|3.13|) 
and ( |3.43| ) it follows that (G + (p)) q can be expressed via the two rational functions T(p) 
and R(p) in the following way 



(G + ( P )) q = t( p ) - ^T( P y - (-iy+*R(py (3.44) 

where the negative sign in front of the square root ensures the correct asymptotic 
behaviour of G + (p) as p — > oo. Now the requirement that G(p) must be analytic in 
the complex plane except from a cut [a, b] implies that the the square root term above 
can have singularities only at a and b and therefore must decompose as T(p)^/p with 
T(p) another rational function. Hence we can parametrize G+(p) in the following way 

G + (p) = = {(P 2 - e 2 ) [A{p)g^{p) + B(p) ] } 1/q (3.45) 

where A(p) and B(p) are polynomials of degree less than or equal to q — 2 and where we 
have made use of the function g_ (p) in order to obtain the property G + (e) = assumed 
earlier. Noting that in the relation (ft.44|) both T(p) and the function appearing under 
the square root, for given I and q, are functions of a definite parity one finds that the 
same must be true for A(p) and B(p). More precisely 

A(-p) = (-l) l+1 A(p) } B(-p) = (-l) l B(p). (3.46) 
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To determine the polynomials A(p) and B(p) as well as the parameter e it suffices to 
evoke the relation fl3.13| ) which implies 

(_!)*+' (p2 _ e2)?-i = (pa _ e 2 )jB 2 (p) _ ^2 _ ^ A 2 (p) _ 2 ^ p A{p)B{p). (3.47) 

and where we note that the number of equations exactly matches the number of un- 
knowns. However, this set of algebraic equations may have many different solutions 
and we must add some boundary condition to select the correct one. Let us note that 
equations ( |3.45|) , (|3.46| ) and ( |3.47|) do not depend on I but only on its parity. We claim 



that the different solutions of equation (p. 47 ) correspond to different values of I. For a 



given I the correct solution can be identified for instance by its asymptotic behaviour 
in the a — > limit. As mentioned in the previous section e(a, b) always lies on the 
positive imaginary axis and in the limit a — > 0, it behaves as a u . More precisely as we 
shall see in section |5] one has 

/ a \ l / q 

e{a,b) ~ lib — (3.48) 



a^o \AbJ 



and this is the criterion which allows us to pick out a unique solution of equation ( |3.47|) . 



One might prefer evaluating the logarithmic derivative p a introduced in (|3.36|) which 
must behave as 

Pa ~" = -- (3-49) 

a^O q 

Let us close this section by considering some explicit examples. In each case the func- 
tion G(p) is determined by the equations ( p^5|) , ( ggBQ , ( $A~7\ ) and |T48|) or flOg ). 



The case 1=1, q=2, i.e. n=0: Here equations (|3.46|) , ( [3.47D and the condition 
that the degree of A and B is less than q — 2 imply that 

a 2 b 2 /— 
A(p) = 1, B(p) = and e 2 = — — i.e. e = +iVab. (3.50) 

e 2 

The expression (|3.45|) for G + (p) hence reads 

G+(p) = = hp - p^] =={Vp~ iPib ~ «)} 1/2 • (3-51) 

After performing the transformation ( |3.3| ) one finds the familiar form of the solution of 
the saddle point equation of the hermitian 1-matrix model with behaviour G{p) — > 
p — > oo 

G{p) = ^ ( 3 - 52 ) 

^(p-a)(p-b) 
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Furthermore one easily verifies that with the expression ( |3.51| ) for G + (p) the formulas 
in section ( |3.2| ) correctly reproduce the usual contour integral representation of the 
solution of the 1-matrix model. 



The case 1=1, q=3, i.e. n=l: Let us emphasize that this set of models contains 
the Ising model on a random lattice as a special case. Since the polynomials A(p) and 
B{p) are of degree less than or equal to q — 2 = 1 and obey the parity condition (|3.46|) 
we write them in the following way: 

A(p) = c, B(p) = p (3.53) 

The constant c and the parameter e are determined by (|3.47|) . For c one finds 

-2e^/e 



e 2 — a 2 b 2 /e 2 

while e is given by 



(3.54) 



e 2 = -eab, e 4 - 6e 2 - As (j + - J - 3 = (3.55) 

According to ( |3. 48| ) and ( |3.49| ) we have to choose the branch of the solution of this 
fourth degree equation for which e > and p a — > 1/3 when a — > 0, i.e. when e — ► oo. 
For the solution which matches these criteria one has 



where it is understood that the positive square root should be taken. 

3.4 The two-loop correlator at genus zero 

One way to calculate the two-loop correlator is to use the directly the recipe 

However, there exists a less work demanding method. The two-loop correlator at genus 
zero must satisfy the following equation 

W {p + iO,q) + W°{p-iO,q)+nW°{-p,q) = - . 1 - , pe[a,b] (3.58) 

\P-Q) 

which appears when one applies the loop insertion operator d/dV(q) to the saddle 
point equation (|2.16|) . This is an equation of the same type as fl2.16|) . One can split 
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the two- loop correlator in a regular and a singular part. The regular part is easily 
found and coincides with what one finds by acting with the loop insertion operator on 
the regular part of the one-loop correlator. The singular part fulfills the homogeneous 
version of the equation ( 3.58[ ). To solve this equation it is convenient to perform a 



rotation like ( |3.2| ) for each of the two variables of W(p, q) so that one has 

W(p, q) = -e^W ++ (p, q) + W + „(p, q) + W- + (p, q) - e~ iun W — (p, q). (3.59) 

with 

W + .(p,q) = W—(-p,q) = W-+(-p,-q) = W ++ (p,-q). (3.60) 

We note that 

Wy(p,9) = -J^ W M)i hJ G {+>-}■ (3-61) 



Now the singular part of W® + (p,q) fulfills an equation similar to ( [3. 11 ) in each of 



the variables and a parametrization of the most general solution can be written down 
using the functions G + (p) and g + (p)G + (p) introduced in section ( |3.2|) . The following 
requirements on W(p, q) single out a unique solution. 

• W(p, q) = d y( p j dv(q) ^ mus ^ be symmetrical in p and q and regular when p = q 

• W°(p, q) can have a singularity of the form ((p — a)(p — b))^ 1 ^ 2 but no additional 
poles at a or b since W°(p) has only a singularity of the type ((p — a)(p — b)) 1 ^ 2 . 

• W(p, q) has the following asymptotic behaviour 

W{p,q) ~ 0(l/p 2 ), p^oo. (3.62) 

The unique solution reads 

pg+(p) - qg+iq) 



W° ++ (p,q) = -±-^\G + (p)G + {q) 



-1 — a- 



(P9+(p) ~ Q9+(q)) [pVp ~ q\fq 



p2 _ q2 



(p 2 — q 



2\2 



(p + q)' 



(3.63) 



We see that the result does not show any explicit dependence of the matrix model 
potential. Hence the universality of the two-loop function observed for the 1-matrix 
model Jl6|, 0, |18| extends to the 0(n) model on a random lattice. In addition there is 
a large degree of universality with respect to n. (We remind the reader that the result 
above is valid for any n, but that different values of n give rise to different functions 
G+(P)0 
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In accordance with the fact that W+(p) depends on the potential V(p) only via 
W r -(p) and that 

w^-JL-'fJ-) (3.64) 



dV+(q) 4 — n 2 dq \p — q 

we find that the two-loop correlator can be written as a total derivative 

WM = (-4- + G + (p)G + (,)75 gg±MzLj|±M \ . (3. 65 ) 

4 — n A oq [ p + q p z — q z J 

To proceed with the solution of the loop equation we need to know W®_(p,p). To 
determine this quantity we must analyse carefully the limit p — > q of W+ + (p, —q) 
(which is a rather time consuming task). The outcome of the analysis is 

n 1 fl e 2 — a 2 1/ 2 i2\2 P 2 1 

W +-^ P) = A^T 2 \2{p 2 - a 2 ){p 2 - 6*) + 4 (a ~ h) (p 2 -a 2 ) 2 (p 2 -b 2 ) 2 ]- {3M) 

We draw the attention of the reader to the fact that W+-(p,p) is a rational even 
function with poles at p — ±a and p = ±6. This will be of importance for the 
following. 

3.5 The one-loop correlator at genus 1 
3.5.1 The structure of the 1-loop correlator 

For genus 1 the loop equation reduces to 

j^W (1 Xp) = W°_(p,p) (3.67) 

where is the linear operator entering the left hand side of the loop equation (|3.6|), 
i.e. 

4^a/(p) = [K + ~ W° + (p)} /_(p) + [K_ - W°(p)} U{p). (3.68) 

Let us note for later convenience that using the decomposition ( p.7| ) we can write the 
action of the operator K on & function f(p) as 

Kf(p) = (4-n 2 )/ ^^L^{ Wa+ (u)f.(u) + Ws-(u)f + (u)} 
Jc 2 2m p A — uo z 

= 2 (4 — n 2 ) ■ even fractional part of W s+ (p)f-(p). (3.69) 



Noticing the structure of expression ( |3.66|) for W9_(p,p) and bearing in mind the 



strategy for calculating higher genera contributions in the case of the hermitian 1- 
matrix model we will seek to express W l (p) in the following way 

W l (p) = 4 2 M 2) (P) + B® X ?\p) + 4 l) xi l \p) + B[ l) xl l \p) (3.70) 
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the idea being that the x-functions should allow us to invert the operator K, i.e. 

^(p) = Ty^p Kxt\p) = (-T^p (3-71) 
We will require that the x-functions have the following asymptotic behaviour 

X { r\p),X ( r\p) ~ o(-\ (3.72) 



p^oo 



p2 



where the possibility of a 1/p term has been excluded in order to ensure that the 
relation ( [2.14| ) remains true. The coefficients and b[ 1 \ i — 1,2 follow from the 
decomposition of W+_(p,p) into fractions of the type (p 2 — a 2 )~~ m and (p 2 — b 2 )~ m , 
m = 1, 2 and read 

4 1} = ^ (e 2 - « 2 - ^(a 2 + 6 2 )) , A? = \a 2 , (3.73) 
S« = (e 2 - a 2 - i(6 2 + a 2 )) , = \b 2 . (3.74) 

3.5.2 Determination of the x-fimctkms 

Since the analyticity structure of the x-f unc ti° ns should be compatible with that of 
the 1-loop correlator it is natural to try to construct these functions using as starting 
point the functions G(p) andpG(p) (cf. to equation ( j3.19|) ). From ( |3.69| ) it follows that 

KG(p) = KpG(p) = 0. (3.75) 

Next, let us consider the following functions 

^M = (^. = (*™> 

Applying the operator K to these functions one finds 

(P 2 -a 2 )' a ^(p 2 -a 2 )' v ; 

where {m^j} and {w^ } are some constants. This means that from either of the two 
series of functions (p) and (p) we can construct functions Xa (p) obeying (|3.71|) 
and ( p.72j ). (We remind the reader that G(p),G(p) ~ 0(l/p), p — > oo.) However, nei- 



ther of the two series alone can serve as building blocks for x^ (f>) since all the functions 
<f>w (p) and (p^ (p) have poles at p = —a which contradicts the assumption concerning 
the analyticity structure of W(p). We are henced forced to take linear combinations 
of 0's and (/>'s to kill these unwanted poles. One type of such linear combinations with 
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correct analyticity properties leaps to the eye; the functions (gfs) G(p) ( c f- to equa- 
tion (|3.35|) ). Let us try to build X^(p) from such functions. The expression for 
appeared in equation ( |3.35| ) and for d d u2\l we find using (|3.35|) and ( |3.40|) 



d V „, , 3 1 f d , „ „, J „ ,9 



where the constants C\ and C2 are given by 

r 2 J 1 n , 19 Pa 1 1 9A a \ 

Co ° = 3l^ (1 -^-2^-2A: (1 -^W/ (3 ' 79) 

c * = hw- (3 - 80) 

We note that we always have a recursive relation like (|3.78|) relating the (k + l) th 
derivative of G{p) to the k th and the (k— 1) . This follows from the fact, already evoked 
several times, that any solution of the saddle point equation Q2.16 ) can be parametrized 
in terms of any two other independent solutions. The nature of the parametrization 
follows from an analysis of the analyticity structure and the asymptotic behaviour of 
the functions involved and a recursive relation for the expansion coefficients can be 
found. Since for the moment we will need only an d we shall not enter 

into a detailed discussion of this point, but we will make use of such considerations in 
section |0] concerning the calculation of higher genera contributions. 

Now, let us consider the action of the operator K on the functions above. One finds 

where the moments M.\ and .M2 are defined by 

= G(w)) i= i j2 . (3.83) 

This means that we can choose our x a -f unctions in the following way 

-G{p) (3.84) 



Mi = (4-n 2 )/ ^ulWr+iu) G4u))+W r 4u;)((^) G + (uj) 





1 


d 




da 2 




1 

Moo 





G{V) ~ T^^G(p) (3.85) 



da 2 / M„o da 2 
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Needless to say that xt^iP) appears from Xa"^ip) by the replacement a <-> b. Now 
combining the relations ( |3.73| ), ( |3.74j ) (|3.84j ) and ( |3.85| ) one has an explicit expression 



for the 1-loop correlator at genus one. 
3.6 The free energy at genus 1 

To determine the free energy at genus one we use again the strategy of reference [jTJ; 
namely we seek to express the basis vectors xl" 1 "* (p) and xt"^ (p) as total derivatives with 
respect to the loop insertion operator. The case m — 1 is relatively simple. Using the 
relation ( |3.64| ) as well as ( |3.31| ), ( |3.35| ), ( |3.39| ) and ( |3.40| ) one finds from the boundary 
equations ( |3.29|) and ( |3.30|) after a lengthy but in principle straightforward calculation 



{1) lrflogo^ (1) _ l dlogb 2 

The case m = 2 is less simple but due to the appearance of the factor M\ in the 
relation (|3.85|) is is obvious that Xa(p) m ust be closely related to cHog A4.\/ dV(p). By 



explicit computation one finds that this quantity can actually be expressed entirely 
in terms of x 2 a {p)-> xi^G ); Xb \P)i fl2 > ^ 2 ^ ■> an d ff? which is a non trivial result. 
Let us briefly comment on the key relations which ensure this property. (We will also 
need these relations for our discussions in section |3.7| .) Acting with the loop insertion 
operator on AA.\ as usual implies performing an explicit differentiation after the matrix 
model coupling constants as well as an implicit differentiation after a 2 and b 2 . The 
explicit differentiation leads to the appearance of the quantity -§^P^G(p) which using 
the relations ( |3.31| ), ( |3.35| ) and (|3.40| ) can be written as 




(3.87) 

The implicit differentiations lead to the appearance of mixed double derivatives of G(p) 
which with the use of ( |3.35| ) and ( |3.40|) can be expressed in the following way 



d d if l 1 de 2 ) d 



db 2 da 2 




0<2] 4sG(p). (3.88) 



e 2 — b 2 da 2 J db 2 
In total one ends up with the following expression for x^ (p) 

2 (2) l dlogjMi 1 d\og\a 2 - b 2 \ l rflogq 2 l dlog|a 2 -e 2 | 

Xa W 2 dV{p) 4 dV{p) 2 dV{p) + 4 dV(p) 1 j 
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Evidently the relevant expression for xl (p) appears from ( p.89| ) by the interchange- 
ment a 2 <-> b 2 . Now inserting the here obtained expressions for the x-f unc ti° ns into 
the expression ( 3.70|) for W l (p) we can do the integration and obtain The result 



reads 



Fi = — -log .Mi - — log Jx - -log la 2 - b 2 \ 
24 24 6 

1 11 1 

+- log \a 2 - e 2 \ + - log a 2 + - log \b 2 - e 2 \ + - log b 2 

+fl (j) (3-90) 
where fi f|J obeys the following differential equation 

"^ = < 3 - 91 > 

and where J\ = M.\(a <-* b). We emphasize that this expression for F\ holds for any 
potential V(M) and any n e] — 2,2[. The first three terms of ( |3.90| ) have a structure 



similar to the terms which appeared in the case of the 1-matrix model and one can 
easily verify that the 1-matrix model (n=0) result is correctly recovered. 

3.7 Higher genera and multi loops 

Having calculated W l {p) we are in a position to further iterate the genus expanded 
version of the loop equation (|3.6|) . While the moments and basis vectors introduced in 
section |3.5.2| certainly lead to simple expressions for the genus one quantities presented 
there, they do not give the optimal parametrization of the model when it comes to the 
representation of higher genera contributions. Let us describe now what we consider as 
the optimal parametrization of the model. We will still work with a set of x-functions 
satisfying the relations (|3.71|) and (|3.72 ). However we will change the set of basis 
functions and moments. 

As basis functions we shall use instead of the functions G(p), 

a set of functions {G^\p), G[ k \p)} defined by 



1. G^\p) and G[ k \p) satisfy the homogeneous saddle point equation (|2.16|) 

G {k \p + iO) + G {k \p - tO) + nG (k \-p) =0, p £ [a, b}. 

2. G {k \p) and G ( b k \p) behave near the end points of the cut [a, b] as 



G^\p) ~ (p - a)- k '^(p - b)^ 2 Gf\p) ~ (p - b)- k -" 2 {p - a)" 1 / 2 . 



20 



3. G <yk \p) and G^\p) are analytical outside the cut (especially near —a and —b). 

4. G^\p) and G < b k \p) have the following asymptotic behaviour 

G?>(p),Gf (p)~^, p^oo. 

Here the conditions 1 and 3 ensure that the analyticity properties of G^'(p) and 
G[ k \p) are compatible with those of the one- loop correlator and that KG^\p) and 
KGf\p) will be even rational functions with poles at p — ±a and p = ±6 respectively. 
The purpose of condition 2 is simply to relate the degree of the poles to the index 
k. The conditions 1-3 are satisfied by many different families of functions, but only 
one family of functions fulfills all four conditions. We note that the set of functions 
( -)p2) G(p), (-J2 ) G{p) > introduced in section ^5| fulfills the conditions 1-3 but not 



condition 4. Furthermore we note that for n = 0we reproduce exactly the basis func- 
tions used in reference 0. Now for n = 2cos(z/7r), let us denote by {^'(p), &f \p)} 
the basis functions corresponding to n = 2 cos ((1 — v)^)-, i.e. 

{GP(p)A%)}„ = {G^( P ),G^(p)} iu (3.92) 
These are the functions that will appear in our definition of moments, namely we define 

M k = 2<f p,V'(u)Gi k \u), J k = M k (a^b). (3.93) 
JCi 2m 

In the case of the 1-matrix model (n=0) we have {& k \p), G { b k \p)} = {G^(p), G { b k \p)} 
and hence we reproduce with our definition exactly the moments (up to a factor 2) 
used for the 1-matrix model 0. Furthermore as in the one- matrix model case, we can 
write one of the boundary condition, (|3.29|) , as 



M = (3.94) 

However, we stress that n = is a very special case. In general we will have G^ k \p) ^ 
& k \p). We draw the attention of the reader to the importance of the condition number 
4 in the definition of the basis functions. With a boundary condition of this type we 
will have for a potential of degree d that M k = J k = for k > d — 1. This gives the 
parametrization of the model in terms of the smallest possible number of moments[]. 

4 One may argue that we still have one parameter too much since for a potential of degree d 
described by d coupling constants, we have d — 1 moment variables plus the two variables a and b. 
Indeed there is a constraint that would allow us to reduce the number of parameters by one, namely 
the fact that the free energy has to be dimensionless. Hence, if we defined our moment variables to 
be dimensionless our results would depend on a and b only via a/b. However, we have not found that 
such a redefinition leads to any simplification from a computational point of view. 
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The functions |(^) G(p), G{p) j behave as 0(l/p 2 ) as p — ► oo. Hence with 

these functions as basis functions we would have the unpleasant situation that even 
with a potential of finite order we would have an infinite number of moment variables. 
We have kept these less pleasant moment variables in section |3.5| and section |3~B" 



since 



they render the formulation of the idea of our iterative procedure more comprehensible 
and since they give a particularly simple representation of the free energy for genus 1. 
The nature of the prefactor of 1 /p k+1 in the requirement on the asymptotic behaviour of 
the basis functions is not important for the argument above. However, it is convenient 
for the analysis of the critical behaviour of the model that this prefactor is independent 
of a and b. We choose it equal to one for simplicity. Let us mention that we can also 
write 

M k = {A-n 2 )l ^{^^(^-^(^(w)} (3.95) 

= (4-n 2 )/ ^{^HGSM-H/^HG^M} (3.96) 

The expression ( |3.96|) is particularly appealing since the integrand does not have any 
cut but only singularities in the form of poles at ±a. This follows from the fact that 
Ws±(p) is a solution of the saddle point equation ( |3.11| ) while G*£±(p) is a solution of 
the same equation with v — > 1 — v and that G^'(p) ~ (p — b)~ l l 2 for p ~ b while 
W s (p) ~ (p — 6) 1 / 2 for p ~ b. Since in addition the integrand is odd, the contour Ci 
can be deformed to a small loop encircling the point a. Similarly, in the case of the 
7^ moments the contour C<i can be deformed into a small loop encircling the point b. 
This observation will prove very useful for our considerations in section [5] concerning 
the scaling limit of the model. 

3.7.1 Recursion relations for G^\p) : G^\p) 

From section 13.31 it follows that 



Furthermore from the defining conditions 1-4 one can conclude 

since dG g a ^ fulfills the requirements 1-3 for G^ +1 ^ (p) and the appropriate asymptotic 
behaviour can be obtained by multiplication by a constant. The expressions for 
and A^ 1 -* can be extracted from the relations (|3.35| ) and ( |3.79| ) respectively. One finds 

A^ 0) = 2tan(z/7r/2) A„, (3.99) 
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We will now derive a set of recursion relations which allow us from the knowledge of 
G<W(p), A<°> and A« to calculate any G^(p), k > 2. First we use the fact that any 
solution of the saddle point equation ([2.1 6|) can be parametrized in terms of any two 
other independent solutions to write 

G^\p) = — L_ {Gt%) + c?>G?>(p)} • (3-101) 

Here the prefactor l/(p 2 — a 2 ) generates the correct leading singularity as well as the 
correct asymptotic behaviour of G^ +l \p). The constant c^> is determined by the 
requirement that G^ +1 \p) should not have a pole at p = —a, i.e. 



-i 



c ( fc ) = (3.102) 



Next, by combining ( |3.98 ) and ( |3.101| ) we obtain the following relations between the 
coefficients c^> and A^ 

ci fe) A« = (k + i), \^ - = (3.103) 

From our knowledge of G^ (p) , A^ -* and A^ 1 ** we can now by means of (|3.101|) and (|3.103|) 
easily write down an explicit expression for any G^(p) (and similarly for G^\p)). 
Furthermore it is obvious that the G-functions appear from the G-functions by the 
substitutions v — > 1 — v and we will use for the relations involving G-functions the 
same notation as above just with all quantities being equipes with a tilde. We have in 
addition the following relation between the G and G functions 

PG? } (P) = Gf\v) + 8<PGP(p) (3.104) 
The argument goes as above and the constant sffi is given by 

(*) _ Gj^jO) 

a ~ Gi k \0) ' (3 ° 5) 

By inserting ( 3.104j) in ( |3.98|) and (|3.101|) we find the following expression for 

si k+1) = 4 fc) 4)- (3-106) 
Hence it suffices to calculate s£> . It reads 
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Now, if we define another set of moment variables by 

M k = 2<f c ^uV'(u)G^(u), J k = M k (a^b) (3.108) 

we have 

M k = M fc _x + s^M k (3.109) 

and inversely 

M k = -^M k - ^^M k -i + ■■■ + w 1 (/ l (3.110) 

To derive this equation one explicitly makes use of the fact that M = 0. These 
two relations allow us to move freely between the two sets of variables. However, we 
stress that it is the M-moments which are the fundamental quantities since these, 
as mentioned earlier, give the parametrization of the model in terms of the smallest 
possible number of moments. Working with the M-moments would for a given potential 
of finite degree (or for a given multi-critical point) lead to the appearance of one 
additional parameter. (For a potential of degree d we will have M k = only for k > d 
while M k =0 for k > d - 1.) 

3.7.2 Recursion relations for X^\p) an( i Xb\p) 

We remind the reader that the aim of introducing the basis functions was to be able 
to invert the operator K. Let us therefore examine the effect of acting with K on such 
a function. One finds 



where fi k ,i and r k> i are defined by 

m = (4-n 2 )/ ^(^ 2 -a 2 )'{^ s+ MG^M + ^_MGSH},(3.112) 

JC 2 27T2 k ' 

r k ,i = Hk,i(a <-> b). (3.113) 



From (|3.111|) we can write down a recursive relation for the x-functions, namely 



Xi fe) (p) = — \g^\ P ) - EM*,i-iX?(p)l (3-114) 
A*fc,fe-i I i=i J 

and similarly for xt (p)- ^he /i-coefficients can be expressed in terms of the moment 
variables and the c^'s. One has 

Hk,i = 0, I > k, n kt o = M k - (3.115) 
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and the remaining /i-coefficients then follow from the recursion relation 



l^k+1,1 = /"fc-u-i + ci fc Vfc,/-i (3.116) 

which is a simple consequence of ( |3.101| ). We note that in particular we have 

k— 1 

Hk,k-i= Rc^Mi (3.117) 

8=1 

3.7.3 The one-loop correlator at genus g 

In analogy with what was the case for the hermitian 1-matrix model we have the 
following representation for the genus g contribution to the 1-loop correlator. 

30-1 

W 9 (p) = E {4 m) xi m) W + B^ X { r\p)} (3.H8) 

m=l 

where the x-functions are given by ([3.1 14|) and where the coefficients take the 
form 

4 m) = £ h ) ^■■■^■■■ J ~, (3 . 119) 

with the indices being restricted by the conditions 

(/-/?) + (* - 7 ) = 2-2g, (3.120) 

s I 

E(A-l)+E(Ti-l) < Sg-m-l. (3.121) 

i=i i=i 



That the equation (|3.118|) holds can be proven by induction using as the starting point 



the expression obtained earlier for W 1 (p). Obviously the proof consists in showing that 
with the representation ( |3.118j ) valid for g' — 1, . . . , g — 1 the right hand side of the loop 
equation ( |3.6|) can be decomposed into fractions of the type (p 2 — a 2 )~ m , (p 2 — b 2 )~ m , 
m = 1, . . . , 3g — 1 with appropriate coefficients. Let us just draw the attention of the 
reader to a few essential ingredients of the proof. 

As regards the first term on the right hand side of the loop equation the existence 
of the above mentioned decomposition follows from the fact that the basis functions 
fulfill the homogeneous saddle point equation and the analyticity requirements 2 and 
3 on page This means that a function of the type (jp)G_ (p) can not have any 
cut but must be a rational fraction with poles at p — ±a and p = ±b of order less than 
or equal to k + m. 

The important step in proving that the second term on the right hand side of (|3.6|) 
indeed takes the desired form consists in showing that dMk/dV(p) and dJk/dV(p) can 
again be expressed in terms of basis functions and M- and J-moments. However, 
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due to the relations ( p. 1U9| ) and ( |3.11U| ) it is equivalent to show that dM k /dV(p) and 
dJk/dV(p) can be expressed in terms of basis functions and moments of the type Mi, 



Ji. For simplicity we shall here take the latter line of action. From the definition ( 3.108 
it follows that 

dM k d / < h \, s \ da? , th ^— db 2 f duo Tr( . .dG^(u) , 

where the first term comes from the explicit differentiation after coupling constants 
and the two others from the implicit differentiation after a 2 and b 2 . Exploiting the 
fact that p k+l G^\p) depends on p only via p/a and p/b we can rewrite the first term 
of (|3~T22D as 

| (p G ?)(p)) = -*G?) W - 2a^-^El _ 2 ^jSM (3 . 123 ) 

Furthermore the analyticity properties of — g fe2 allow us to conclude that we have a 
decomposition of the following type 

= ^Gl'V) + E<iG?»(p) (3.124) 

where the f^'s are some constants. From ( |3.88| ) it follows that for k = 1 we have 

« = -S4felog(^41 (3.125) 



2\^da 2 b \b 2 



a 



and the remaining f-coefficients can be found by repeatedly use of the k = 1 relation 
and the relation ( p.98| ). In conclusion one can write dM k /dV(p) as 

d~M k . M fc+ i 



-2*G« - Aa 2 X^ G^\p) - ^G^\p) 



dV(p) a a { a ^ J M 1 

-^E^j^W-^Gi 1 )^)} (3.126) 

where we note that the vf \ terms have cancelled. Collecting the here given information 
it is straightforward to complete the proof of the representation (|3.118|) for W 9 (p). 
In case of the ordinary one-matrix model one has fp'™. p 7 («, &) = ( a ~ b)~ s where 
5 = Ag — 2 — m — J2t=i(Pi ~ 1) _ Z)j=i(7j — !)• In the general case this is no longer 
true. However, we emphasize that we still have that all explicit dependence on the 
matrix model coupling constants is hidden in the moment variables. The function 
fp^ 1j i( a i7 ( a ) b) is a function of the endpoints of the cut only and expressed in terms of 
the variables e and a it takes the same form for all values of n e] — 2, 2[. Unfortunately 
we have not been able to write down the generic expression for ^(e, a). 
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3.7.4 Multi-loop correlators 

From W g (p) we can obtain W g (pi,p2, ■ ■ ■ ,p s ) for any s by repeatedly use of the loop 
insertion operator (cf. to equation Q2.6| )). Analyzing the structure of the loop insertion 
operator, one can write down formulas similar to ([3.118 ) for the multi-loop correlators. 



We will not pursue this aim, but let us mention that from the discussion in the previous 
section it follow that the genus g contribution to the s-loop correlator as in the 1-matrix 
model case depends on at most 2 x (3g — 2 + s) moments for g > 1.. The same statement 
is true for g = provided s > 3. This can be seen from the expression ( |3.57| ) for the 
two-loop correlator at genus zero. We note that the expression ( |3.57| ) could also have 
been obtained by applying the loop insertion operator to the one-loop correlator at 
genus zero. However, this method of calculation is more time consuming than the one 
actually used. 

3.7.5 The free energy 

From W 9 (p) we can obtain F g by application of the inverse loop insertion operator, 
the inversion being possible due to the relation (p. 126|) . One easily infers that as in 
the 1-matrix model case the genus g contribution to the free energy for g > 1 depends 
on at most 2 x (3g — 2) moments and that for g > 2, F g will be a sum of terms of the 
same type as those entering the relation ( |3.119 ) where the indices fulfill (|3.120|) as well 



as a relation like ( |3.121 ) where on the right hand side 3g — m — 1 is replaced by 3g — 3. 



4 The cases n = ±2 

The cases n = ±2 pose no particular problems. On the contrary they are in a certain 
sense easier to solve than the generic cases, namely the saddle point equation as well 
as the loop equations can be expressed in terms of functions of a definite parity and 
the generic solution to the saddle point equation can be parametrized using only one 
singular function. 

4.1 n = -2 



Let us start by noting that for n = —2, if we introduce Aj = in the integral ( 2.8 
we find 

N 

Zoc fn^^WrKA.-A,.) 2 (4.1) 

Jo i=i i<j 

Hence the partition function looks very similar to the one of the usual hermitian 1- 
matrix model. There are two important differences though. Firstly the interval of 
integration is restricted to the positive real axis. While this does not give rise to any 
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complications concerning the solution procedure it shows that the present model clearly 
contains other critical points than the usual hermitian 1-matrix model; namely points 
for which the eigenvalue distribution exactly touches the origin. In this respect the 
model is very similar to the complex matrix model which is given by an integral of the 
same type, the Aj's playing the role of the positive eigenvalues of a matrix 0^0 fl9|j . 
However, there is an important feature which differentiates the 0{— 2) model from 
both the complex and the hermitian one matrix model. The potential V(y/Ai) might 
contain half integer powers of Aj. Likewise the correlation functions that one would be 
interested in calculating will typical involve half integer powers of Aj. Let us proceed 
to discussing how the usual iterative procedure can be adjusted to these circumstances. 

4.1.1 The one-loop correlator at genus zero 

As in the previous sections we will assume that the 1-loop correlator W(p) (defined 
by ( |2.6| )) is analytic in the complex plane and that it behaves as 1/p as p — > oo. Let 
us decompose W{p) as 

W(p) = W+(p) + pW^(p) (4.2) 

where the functions W + (p) and W-(p) are both even in p. Now W + (p) and W-(p) have 
in addition to the cut [a, b] a cut [—b, —a] and the analyticity requirement on W(p) 
implies 

W + (p + iO) - W + {p-iO) =p(W„(p + iQ) -W_(p-iO)), pe[a,b) (4.3) 

In particular the eigenvalue density can be found from either one of the two functions 
W + (p) and W-(p) (cf. to equation (|2.15|)) 



p(p) = JL (w_ (p + iO) - W- (p - iO)) , p G [a, b] 

VK 

= — (W+(p + iO) - W+(v-iO)), pe[a,b] (4.4) 
in 

The saddle point equation ( |2.16|) becomes an equation for W- (p) and expressed in terms 
of the variable p 2 instead of p it takes the same form as the saddle point equation of 
the hermitian 1-matrix model. Hence the solution of the present equation can be read 
off from the solution of the latter. One finds [Q 

w w = i / ^n-L ( (r 2 - ° 2 >(? 2 - " 2 > T ,4.5) 

whew a 2 and b 2 are given by 

- °> t 7T~T-,. ~F>T, .9 uKTF). = 2 ( 4 - 6 ) 



T Cl 2ni(uj 2 -a 2 y/ 2 (co 2 -b 2 ) 1 / 2 ' J Cl 2m (oo 2 - a 2 y/ 2 (u 2 - b 2 ) 1 / 2 
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We note that from W-(p) we can find W + (p) by the following recipe 

W + ( P ) = 2 f *tL W+ (u>)-^ = 2 I ^Lw.^^s (4.7) 
JCi 2m p A — uj z JCi 2m p z — 

4.1.2 Higer genera and multi loops 

Let us introduce a decomposition of the loop insertion operator, namely 

d d d . . 

+ P „ r/ \ (4-8) 



dV(p) dV+(p) ^dV_(p) 

where the operators d/dV + (p) and d/dV_(p) contain only even powers of p. Then we 
can rewrite the loop equation ( |2.22j ) as 

JL i tovps = f + i i w _ M (4 . 9) 

p z Jci 2m p l — u A M z dV_(p) 



where we have explicitly made use of the relation (|4.3| ). Instead of searching a solution 
of the equation (|4.9|) one can search a solution of the following equation 

L s3 B '- M = {w - wf + w^ w - { '" ] (4 - 10) 

since such a function will automatically fulfill 

duV(u)W-(u) = (4.11) 

Ci 

The genus g contribution to the free energy of the 0(—2) model now takes the same 
form as the genus g contribution of the free energy of the hermitian one-matrix model 
given in reference M provided the moments and Jk are defined by 

r duj V'{uj) 
Mk ~ J Cl 2m{uo 2 - a?) k +y 2 {uo 2 -b 2 ) 1 ' 2 ^ 
J k = M k {a 2 ^b 2 ) (4.13) 

and the parameter d is replaced by 

d = b 2 -a 2 . (4.14) 

This statement is easily proven. First one rewrites the loop insertion operator d/dV-(jp) 
in the moment parametrization and realizes that it takes the same form as the loop 
insertion operator of the hermitian one-matrix model (with the modifications given 
above) except for p being replaced by p 2 . This means that the analogy between the 
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loop equations of the two models holds to all orders in the genus expansion. Secondly 
one notes that for the 0(— 2) model the following obvious relation holds 

mM = im F ' (4 - 15) 

and the correctness of the statement concerning the free energy becomes evident after 
a few moments thoughts. 

We emphasize that the contour C\ above encircles only the cut [a, b\. If the potential 
is even, however, one can immediately rewrite the integrals above as integrals along the 
contour C<i- Then performing the change of variable uj 2 —>■ uj one reproduces exactly 
the expression for the free energy of the hermitian one matrix model (of course with 
the assumption that the support of the eigenvalue distribution lies on the positive real 
axis) . 

We will not pursue the explicit calculation of multi-loop correlators for the 0(— 2) 
model in the present publication but let us emphasize that such calculations pose no 
particular difficulties. One simply rewrites the loop insertion operator in the moment 
parametrization, using the boundary equations ( |4.6j ) and applies it to the free energy. 
As mentioned above d/dV_(p) has a structure similar to the loop insertion operator of 
the hermitian 1-matrix model. The even part d/dV + (p), however, is less simple and 
involves elliptic integrals. 

4.2 n = +2 

4.2.1 The one-loop correlator at genus zero 

Let us introduce again the decomposition of the 1-loop correlator given in equa- 
tion ( |4.2j ). As before we then have the relation ( |4.3|) between W + (p) and W_{p) and as 
before the eigenvalue density can be found from either of the two as described in equa- 
tion (|4.4j). The saddle point equation turns into an equation for W+ip). This equation 
when expressed in terms of p 2 takes the same form as the saddle point equation for 
the hermitian 1-matrix model and the solution of the present equation can be found 
by exploiting the analogy with the latter. The result for W+(p) reads [|J 

1 r (tf-aW-Vn* 
+ KP > 2Jc 1 27tip 2 -u 2 \{uj 2 -a 2 ){uj 2 -b 2 )f y ' 

Of the two boundary conditions which determine a 2 and b 2 one ensures the correct 
asymptotic behaviour, W + (p) ~ 0{l/p 2 ) as p — > oo, and can be written in the standard 
form 

du V'{uj)uj 



d 2m (uj 2 -a 2 y/ 2 (uj 2 -b 2 ) 1 / 2 



(4.17) 
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The other one expresses the fact that the eigenvalue distribution is normalized to one 
and reads 

iM w >^\ (418) 

or 

- t dp(v 2 -a 2 ) l l 2 (b 2 -p 2 ) 1 / 2 I —— ^ = 2 (4 19) 

As opposed to what is normally the case this second condition can not be written as 
a single contour integral. This is due to the fact that W+(p) contains only the even 
powers of p, i.e. the behaviour W(p) ~ 1/p can not as usual be imposed by simply 
referring to the contour integral ( |4.16| ). 

Even though the complexity of the second boundary equation does render the iter- 
ative calculation of the free energy and the multi-loop correlators more involved than 
for n = —2, the moment technique is still applicable. However, a detailed analysis 
of the structure of the free energy and the multi-loop correlators at higher genera is 
rather work demanding and we shall in the present publication restrict ourselves to 
exemplifying the applicability of the moment description by calculating the free energy 
at genus 1. Our line of action will follow closely the one taken for n g] — 2, 2[. 

4.2.2 The two-loop correlator at genus zero 



Introducing the decomposition ( |4.8| ) of the loop insertion operator we can write the 
loop equation ( |2.22| ) as 

f ^^P4w + (P) = (W + ( P )f + ^rfinS^O.). (4.20) 
JCi 2ixip A — or iV z dv + [p) 

which in its genus expanded version reads 

[K - 2W° + (p)} WUP) = E w((p)Wl a \p) + ^J—Wf 1 ® (4.21) 

where 

K!(P) = i (4-22) 

JCi 2,7X1 p z — UJ Z 

To proceed with the solution we need to calculate the following two-loop correlator 



The simplest way to do this is to proceed as in section |3.4| . From (|2.16| ) it follows that 
W+ + (p,q) must fulfill the following saddle point equation 

Wl + {p + zO, q) + Wl + (p - iO, q) = -- P ' p G [a, b] (4.24) 
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The complete solution of this equation, consisting of the sum of a particular solution 
and the complete solution to the corresponding homogeneous equation, is easily written 
down. Using the fact that W+ + (p,q) must be symmetric in p and q, finite for p = 
q, have the asymptotic behaviour W9,(p, q) ~ 0(1 /p 2 ) as p — > oo and behave as 
((p — a)(p — b)) 1 ^ 2 in the vicinity of a and b. one finds that it must necessarily take 
the form 

W° ++ (p, ,) = i=L=L L + (f + , 2 ) (l - (^f^) 1 1 (4.25) 




where C is some yet undetermined constant. Now the boundary equation (|4.18| ) implies 
that W ++ (p, q) must fulfill the following equation 

I dpW ++ (p,q) = 0, Vg (4.26) 

JCx 

From this equation one can extract the value of C. This is most easily done evaluating 
the integral at q — 0. The result for C reads 

where 

f h 2_ 2\l/2 

fc » = ( J > ^ = ^ ( 4 - 28 ) 

and where K(k a ) and E(k a ) are the complete elliptic integrals of the first and the 
second kind respectively. To determine W® + (p,p) which is the quantity which enters 
the loop equation we must analyze carefully the limit p — * q of the expression ( |4.25|) . 
One finds 

w o (p p) = I 2 1 P y-& 2 ) 2 

++^'^ 8(p 2 -a 2 )(p 2 -6 2 ) 16 (p 2 - a 2 ) 2 (p 2 - b 2 ) 2 y ' } 

We note that the right hand side of the loop equation ( |4 . 2 1|) for g = 1 takes the same 
form as in the case n g] — 2,2[ , the constant C playing the role of e 2 — a 2 (cf. to 
equation (|3.66|)). 



4.2.3 The one-loop correlator at genus one 

We shall try to express W\(p) as in equation ( [3.70D with the function XaHp) an d Xb (p) 



obeying again the relations ( |3.71| ) and ( |3.8| ) with K given by ( |4.22| ). The corresponding 
A and B coefficients read 

4 X) = \~o~~~U> (C-I(a 2 + b 2 )) , 4 2) = ^a 2 , (4.30) 



a 



b 2 V 2 V V' 16 
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In analogy with the case n e] — 2, 2[ we will express the x-fu nc ti° ns m terms of a set 
of basis functions {G^\ G^}. To begin with let us introduce 

0(O)(P) = (p2_ a 2 ) l/2 (p 2_ 6 2 ) l/2 ( 4 - 32 ) 

This function clearly fulfill the following identity 

K^°\p) = 0. (4.33) 
We now define G^ (p) and Gj® (p) for > 1 by the following requirements 

1. G^\p) and G^\p) are even in p and fulfill the homogeneous saddle point equa- 
tion 

G(p + iO) + G(p - iO) = 0, pe[a,b). (4.34) 

2. In the vicinity of the endpoints of the cuts the functions G^\p) and G^\p) 
behave as 

G^\p) ~ (p 2 - a 2 )-^ 1 / 2 ^ 2 - b 2 )- l '\ Gi k \p) ~ (p 2 - a 2 )- 1 / 2 ^ 2 - ft 2 )-^. 

(4.35) 

3. G^(p) and G^\p) are analytic everywhere else 

4. They fulfill the conditions 

/ dpG { a k) {p) = <f dpGf\p) = 0. (4.36) 

5. They have the following asymptotic behaviour 



G {k \p), Gf (p) = const • 0<°>(p) + -L_ + O , 



p -> oo. (4.37) 



The role of the three first requirements is the same as in the generic case. Condition 
number 4 ensures that the eigenvalue distribution stays normalized to all orders in 
the genus expansion (cf. to equation ( |4.18| )). In the generic case this could simply be 
taken care of by demanding that W 9 {p) ~ 0(l/p 2 ), p — > oo for g > 1, i.e. by excluding 
the possibility of terms of order 1/p in the G-functions. However, in the present case 
we are calculating only the even part of W(p) so the normalization condition must be 
imposed in a different way. Condition 5 is chosen with the aim of rendering the moment 
variables as simple as possible. This should become clear shortly. The conditions 1-5 
determine the G-functions uniquely. One has 

G {k \p) = €\p) + S^^\p) (4.38) 
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where ^ 

^\P) = ( p 2 _ a 2)fc+l/2( p 2 _ 62)1/2 ( 439 ) 

and 

and similarly for G^\p). For our considerations in the following section we will need 
the explicit expressions for and They read 

Si" = -Ajl-!^L (4.4D 



a 2 - b 2 1 AT(fc, 



From the basis functions it is straightforward to construct the x-f unc ti° ns - F° r that 
purpose let us consider the action of the operator K on the G-functions. One finds 

k—l i k—1 -i 

KGi k \p) = g M lw —^ v KGtXv) = E Jij^Tpy (4-43) 
where the moments Mi and J; are given by 



M k = £^ MGPM .£^ ( ^_ a2) ^ 2 _ y)i/2 (,44) 
J fc = M fc (a 2 ^6 2 ) (4.45) 



The advantage of imposing the requirement 5 on the G-functions should be clear by 
now. One could have taken G^'(p) as a linear combination of 0w(p) with any <j${jp) 
with / < fc. However, due to the condition ( 4.17Q we obtain a particularly simple 



expression for the moments by choosing 1 = 0. From ( |4.43| ) it follows that the \- 
functions are given by 

Xl k) (p) = ^ |g?>(p) - Y;M k „ l+1 GU(p)\ , xi k \p) = Xi k) (p) (a «-> b) (4.46) 

We note that for a potential of degree p one has M q = J q = for q > p. Now all the 
elements in the representation ( [3.70 ) of the 1-loop correlator at genus one have been 



determined and it is easy, collecting the results of the present section, to write down a 
completely explicit expression for W\{p). 

4.2.4 The free energy at genus one 

To determine the free energy at genus one we use the usual strategy of expressing the x~ 
functions as total derivatives with respect to the loop insertion operator d/dV + (p). The 
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key point in this procedure consists in determining da 2 / dV+(p) and db 2 /dV+(p). These 
quantities can as usual be extracted from the boundary conditions ( [4.1 7|) and ( |2.14| ). 
The actual calculation is more involved than usual but after the use of various relations 
between elliptic integrals one arrives at the following pleasant result 

x <u fp) = dl °z a2 X W fp) = ^! (4 47) 

Having obtained the expressions for da 2 / dV + (p) and db 2 / dV + (p) it is relatively straight- 
forward to show that 

,2 v (2) fnU 2 d log Mi l d\og(a 2 -b 2 ) 2 dloga 2 
Xa U ) " 3 dV+(p) 3 dV+ip) 3 dV+(p) { ' j 

and similarly for \b \Pj- Now combining the A and B coefficients given in ( |4.30| ) 
and (|4.31|) with the here obtained expressions for the x-functions one finds that W^_(p) 
indeed takes the form of a total derivative. The free energy at genus one can hence be 
extracted and reads 

Fi = -^logMi-^logJx-^log^-a 2 ) 

+z^ 2 + h logb2 ~ \ log {K{ka)) ~ i log {K{kb)) (449) 

It is interesting to note the similarity of ( |4.49| ) with the expression (|3.90 ) obtained for 
n e] - 2,2[. 



5 The critical regime 
5.1 The critical points 

As mentioned earlier the matrix integral defining the 0(n) model ceases to exist when 
the support of the eigenvalue distribution approaches zero, i.e. when a — > 0. This 
gives rise to a new set of critical points for which no analogues exist for the 1-matrix 
model H, [H H . These are the critical points that we will consider in the following. We 
will take v to be in the interval < v < 1. Then we airways have a |e| <C b (since 
as we shall see very soon e ~ a u ) which simplifies the analysis. Although not more 
complicated the cases v — 0, 1 require special treatment. 

At the singular points the eigenvalue distribution vanishes at one endpoint of its 
support (here a = 0) with a critical exponent, /i, or equivalently 

W s {p)~p^ p^O. (5.1) 

Let us recall the possible values of \i for the 0{n) model 0, |5|]. These can be read of 
from the expression ( |3.23j ) for W® + (p). Obviously the possibility of new types of critical 
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behaviour is due to the presence of the function G(p). For a = (or equivalently 
a <C p ~ b) the function G(p) takes the form 

G(p) ~ — - cosh-?/;— — — — — — , cosh-?/; = — -. (5.2) 

b sinh^ p 



This is most easily seen by verifying that the the function (|5.2|) satisfies the criteria 
1-3 on page |8|. Now letting p — > we find that 

G{p) ~ p" -1 + ©(p - "* 1 ) and hence G(p) ~ p - " + 0(p"). (5.3) 

It then follows from (|3.23|) , (|3.20|) and (|3.3|) that by fine tuning the potential of our 
model (i.e. the polynomials A(p 2 ) and B(p 2 )) we can reach for a given value of v (or 
n) the following two series of critical points 

/i 2m+1 = 2m + l-z/: A(p 2 )~ P 2m , B(p 2 )~(D( P 2m ) (5.4) 

fi 2m+2 = 2m + l + u: B(p 2 )~p 2m , A(p 2 ) ~ C?(p 2m+2 ) (5.5) 

The possible values of /i are exactly those for which n = — 2cos(/i7r). Furthermore it 
can be shown that 7 s4r = — 2z//(/i + 1 + ^)f| d, ||. When ^ = -, with < Z < q and 
l,q E Z, the critical points being characterized by the exponents \x = 2m + l±i/ exhibit 
the scaling behaviour characteristic of 2D gravity interacting with rational conformal 
matter fields of the type (q, (2m + l)q ± Z). However, the continuum theories that one 
obtains from the 0(n) model do not contain all the operators of the corresponding 
minimial models [|, [5|. 

For later book-keeping purposes, let us arrange all critical points into one series 
where the M'th multi-critical point is characterized by 

Hm = M - r}M+i, mk = v, ??2fc+i = l-^ (5.6) 

We note that this definition reproduces the usual notion of a M'th critical point of the 
1-matrix model [y = |) case. 

5.2 Scaling at a M'th critical point 

In this section we will calculate the scaling behaviour of the basic elements of our 
description, i.e. the functions G^l(p) and the moments {M k , J k }. Knowing the scaling 
properties of these objects we can easily extract continuum results from our exact 
results or develop a procedure for calculating directly continuum quantities. 

The most fundamental quantity of our description is the function G(p). From G{p) 
all other quantities can be derived. One can show that in the scaling region (a ~ p < 6) 

lie sinh vcj) p 

G{p) ~ r . r 7 > cosh^ = — 5 - 7 

ab sinh <p a 
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The prefactor comes from the relation ( |3.13| ). Now matching the expressions 
and ( |5.7| ) in the intermediate region a <^ p <^ b one can determine e to leading order 
in a. The result reads 



V 



2ib ( ) ' A* = V (5i 



which we note justifies our statements concerning e made in section |3.3| . Hence for 
p ~ a we have in accordance with the analysis of the previous section 

G{p) ~ a" -1 , G{p) ~ a"" (5.9) 

Knowing the scaling of e we can furthermore determine the scaling of all G^'(p) and 
Gf\p). Namely, from the relations ( |3799| ) and (|3TT00| ) we see that ~ a 2v ~ 2 and 
A^ 1 ** ~ a~ 2u and then the recursion relations (|3.103|) tells us that 

Ai fc) ~ -4 ~ a ~ 2%+i ( 5 - 10 ) 

In particular G^(p) ~ a~ 1- ^ and in general 

G< fc) (p) ~ a"*-"^ 1 , Gi fc) ~ a - *-"*, p~ a (5.11) 

while all the G^{p)- and Crj; -functions for p ~ a become proportional to G(p) and 
G(p) respectively. 

Let us now examine the scaling properties of the moment variables. We remind the 
reader of the fact that the integrals defining and Jk when written in the form ( |3.96| ) 
reduce to local integrations around a and b respectively. Since for a <C p ~ b, G^ k \p) 
as well as W s {p) are independent of a we have 

J k ~ O(a ). (5.12) 

The M-moments, on the contrary, have a non trivial scaling. By definition of a M'th 
multi-critical point one has at such a point 

W» M {p) ~ p M -w+i _ a M-^ M+w ^ M ^ _ ( 513 ) 

It now follows that 

M k ~ fl tf -*+(w-*) (5.14) 

We see that for < M, scales with a positive power of a and that Mm ~ o°- The 
moments M*. with k > M are equal to zero. This can be seen by deforming the contour 
to infinity. 

Having determined the scaling properties of basis functions and moments it is easy 
to pass to the continuum limit. For instance to determine the genus one contribution 
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to the 1-loop correlator in the scaling limit it suffices to note that in this limit the right 
hand side of the loop equations for W l (p) reduces to 

~ \ / \ 1 (a 2 1 2v 2 - Av + 1 1 1 
KW {P) = 4^ |T(^^ 4 ■ (5 ' 15) 

This immediately tells us that 

W\ P ) = £ X ?>(p) ~ ^ ~1 U + \ ^\P). (5.16) 

We note that there is no simplification of the x-functions in the scaling limit. All 
terms in the relation ( [3.114j ) are of the same order in a. Now using the relations (|3.84|) 



and fl3.85| ) bearing in mind that Ai\ = X^s^'Mi we find the following expression for 



Fi in the scaling limit. 

F 1 = -llog(M 1 a-^ 2+13 ^ 5 ). (5.17) 

The exponent of a vanishes if and only if v = ~. Hence we reproduce correctly the 
1-matrix model result and we see once again that this case is very particular. From 
the expression ( |5.16| ) of the scaling relevant part of W 1 (p) one can pursue the iterative 
solution of the loop equation directly in the continuum. This only requires that one 
writes down a continuum version of the loop insertion operator. Let us mention a few 
properties of this operator. First of all one finds that the loop insertion operator in the 
scaling limit reduces to a differentiation after a 2 and the moments M\~. This implies 
that, not surprisingly, no J-moments will appear in the scaling limit. Furthermore the 
dimension of the loop insertion operator can easily be extracted. It equals a~ MM ~ 2 . 
Let us stress that the expressions ( |5.16| ) and ( |5.17| ) as well as all results that one 
would obtain by further iterations of the loop equations are valid in the vicinity of any 
M'th critical point and independent of which detailed prescription one might choose 
for approaching such a point. However, whenever needed one can easily specialize to a 



given scaling prescription. In section |5.4| we will show how one can calculate explicitly 
the moments when one approaches the critical point by tuning an overall coupling 
constant of the potential. 

5.3 The basis functions in the continuum 

As explained in the previous section the G^-functions do not play any role in the 
scaling limit. Let us write the G^-functions in this limit as 

1 



G { a\p) = ~ 7 l^Tx ^fM (5.1* 

2 cos [vis 1 2) b K+L 
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where we use again the parametrization p = —a cosh <fi. Now fk{4>) is a dimensionless 
function and all dependence on the scaling parameter a is hidden in the prefactor e Uk 
where 

(5.19) 

The value of the exponent Vk follows from the relation (|5.11 ) and reads 

Vk = —k - rjk+i. (5.20) 
The expression for the function /o(0) can be read off from (|5.7|) . One has 

(5 - 21) 

The remaining /^-functions can be found from the continuum versions of the recursion 
relations flpp , (|3TToT|) and ( glQg ). They read 

J fc /*+i(0 = ( 5 - 22 ) 

= is ^2 7 (A-iW+7fc/fc(0)) (5-23) 
lo smh 

and 

/ fc7fc = 32(A; + l/2), lk+i - h-x = ZlMk+i- (5.24) 

The new dimensionless parameters 7& and 4 are related to the original ones and 
A« by 

7fc = lim -i^e" 2 ^ 1 , l k = lim 326Ai fc) e 2% + 1 . (5.25) 

a— >0 o a— >0 

From the expressions ( |3.99|) and (|3.100|) we can determine lo an d h and this enables 
us to solve exactly the recursion relations (]5.24|) . We find 

, A , {Vk*\ Vk(l + Vk)(2 + Vk) ■ ■ • (fc + y k ) f . 

k = (1-V«)(2-V k )...(k-VK) ■ (5 ' 26) 

One can derive additional interesting properties of the quantities appearing above. For 
instance one has 

kk+i = ^§v k {v k - 1) (5.27) 
and it appears that the /-functions satisfy the following differential equation 

/£(«/>) + 2(k + if^f'M + + I) 2 - Vl) fM = 0. (5.28) 

As usual the relevant expressions for the G-functions appear from those of the G- 
functions by the substitution v — > 1 — v and we will use for the relations involving 
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G functions the same notation as above just with all quantities being equiped with a 
tilde. Equation fl3.104j ) relating G- and G-functions translates to the scaling limit as 

-4cosh0/ fc (0) = cot {^pj {/ fc _i(0) +a k f k {(j>)} (5.29) 

where the dimensionless parameter a k is related to s%> by 

CT fe = lime- 2r? *-4 fe) - (5.30) 



Using the relations ( |3.106| ), ( |3.107| ), ( |5.25| ) and Q5.26 ) one can determine a k explicitly. 
It is given by 

<?k{k + ilk) = k- (5.31) 

5.4 Explicit calculations at a M'th multi critical point 

In this section we specialize to a particular prescription for approaching a M'th multi- 
critical point. We replace the potential V(p) of our model by ^4j^ where V c (p) is a 
critical potential corresponding to the critical point in question and where T plays the 
role of the cosmological constant (or the temperature). We now approach the critical 
point by letting T — > T c = 1 and define a renormalized cosmological constant A# 

by ill 

T-T c = a^ I+1 ~ u A R . (5.32) 

where /zm = M — t)m+i- That the power of a appearing above is indeed what is needed 
to make dimensionless can be seen by expanding W(p) around W c (p), considering 
p ~ a and using that 

d(TW(p)) 

Of = G (p). (5.33) 

This relation follows from the fact that the expression on the left hand side fulfills the 
conditions that determined uniquely the function G (p). Now, using the relation ( |5.32| ) 
it follows from (|5.14j ) and (|5.17|) that 




6z/ 2 -15z/ + 7] 

1 — > log Ar. 5.34 

Hm + 1 - v J 

This case is particularly simple. Due to the logarithm we do not need to know the 
explicit expressions for the moments in the scaling limit. However, to determine the 
continuum version of any other quantity such expressions are needed. We shall now 
proceed to deriving these. Our starting point will be the relation ( |3.96| ). We remind the 



reader that the contour integral appearing in this relation reduces to a local integration 
around the point a. Hence we only need to know the integrands in the scaling limit. The 
relevant expressions for the G functions appear from the previous section. However, 
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we shall not make use of their explicit form. It suffices to know that they fulfill the 
following relation with x = p/a = — cosh0 

hfk( x ) - x f k {x) = hfk+i(x). (5.35) 

The scaling limit of the function W s (p) has been determined explicitly in reference ||. 
In the vicinity of a M'th multi-critical point one has 

W B (p) = ~e M -*»+ 1 F M (x) (5.36) 



where 



F M (x) = const ■ x M -™+i (b + f da , ) (5.37) 

w \ 7-oo (cosha) A/ "^/y v ' 

The constant B can be completely explicited as a B Euler function but its precise 

form will not be of importance for the following. The prefactor is non universal and 

depends on the critical potential chosen. From the explicit expression for Fm{x) one 

easily verifies that the following relation holds 

(M - r) M +i)F M (x) - xF' M (x) = const ■ f (x) (5.38) 

Now our moments take the form 



M k = 4 bk+1 M k (5.39) 



with 

dx 



k = (4 - n 2 ) j {F M+ (x)f k+ (x) - F M _(x)fV(x)} = (F M , f k ) (5.40) 

where the contour encircles the point x — 1. From the relations ( |5.38|) and fl5.35|) it 
follows that 

~l k M k+1 — (M — k + t)m — Vk)M k - const ■ (/„, f k ) (5.41) 

and since (fo, f k ) cx 5k,o we see that in accordance with the analysis of section |57| the 
moments M k with k > M will vanish. For 1 < k < M we have 

M k+ i = M - k + 7] M - 7] k 
M k Ik 

Hence we can express all our moments in terms of only one, say M\. This allows us 
to determine any continuum quantity up to a non-universal constant. For instance we 
find for Wi(<j>) by means of flBTTp , flO|), QOlD and ( ^42|) 

W\<P) = const ■ {hh{4>) - (jiM + 1 + 6z/ 2 - 12!/ + 3)A(0)} . (5.43) 



Using the relations ( |5.22|) and (|5.28|) one can easily verify that this results agrees with 
the one obtained in the unitary case within the framework of strings with discrete 



target spaces by S. Higuchi and I.K. Kostov |f20 |. 
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6 Conclusion and outlook 



One interesting conclusion which can be drawn from the obtained exact solution of the 
0(n) on a random lattice is that the model exhibits the same kind of universality with 
respect to the potential as the hermitian 1-matrix model. One must expect this kind of 
universality to occur also for two- and multi-matrix models and the present work can 
be taken as an indicator of how one could make use of this universality in the solution 
of these more complicated models. 

It is also interesting to note that our solution provides an exact solution of the Ising 
model on a random lattice. This gives the possibility of studying spin excitations of this 
model away from criticality. Unfortunately the representation of the Ising model on 
a random surface that one obtains from the 0(n) model has vanishing magnetic field. 
However, it is possible to include a magnetic field by adding a 1/M term to the action 
appearing in equation fl2.1|) . In analogy with this one would expect that in general 
the addition of terms with negative powers of M would enlarge the operator content 
of the continuum theories obtained from the model. It would hence be interesting to 
generalize the moment technique to this situation. 

As mentioned in the introduction our d.s.l. relevant moment description of the 0(n) 
model should allow us by comparison with the corresponding moment description of 
the generalized Kontsevich models to determine which is the precise relation between 
the continuum partition function of the 0(n) model for n rational and the r-functions 
of the generalized kdV hierarchies. We have not completed this analysis but let us 
mention a few observations. First of all we see that for the 0(n) model on a random 
lattice we have in the double scaling limit two series of moments with different scaling 
properties. In general for a r-function of the kdV p hierarchy describing the interaction 
of 2D gravity with matter fields of the type (p,pm — 1), ... , (p,pm — (p — 1)) there 
will appear (p — 1) series of moments with different scaling properties [1I1L [T5| . Hence 
the only models for which we could hope for an exact equivalence are the models 
(p, q) = (3, 3m, — 1), (3, 3m — 2). However, as the example with the Ising model clearly 
shows, not even in this case will the equivalence be exact. 

Another interesting aspect concerning the double scaling limit is the interpretation 
of the continuum theories corresponding to non-rational values of v. For instance, one 
might wonder what the topological interpretation of these models is and if there exist 
integrable hierarchies describing them. 

Finally one can remark that the results that we have obtained are actually analyt- 
ical in v. This might open the possibility of attributing a meaning to the model for 
n > 2 and maybe approaching the question of interaction of 2D gravity with matter 
fields with c > 1. 
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